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Abstract 

This paper presents results of wind-tunnel tests that 
demonstrate a novel drag reduction technique for blunt- 
based vehicles. For these tests, the forebody roughness 
of a blunt-based model was modified using 
micromachined surface overlays. As forebody 
roughness increases, boundary layer at the model att 
thickens and reduces the shearing effect of external flow 
on the separated flow' behind the base region, resulting 
in reduced base drag. For vehicle configurations with 
large base drag, existing data predict that a small 
increment in forebody friction drag will result in a 
relatively large decrease in base drag. If the added 
increment in forebody skin drag is optimized with 
respect to base drag, reducing the total drag of the 
configuration is possible. The wind-tunnel tests results 
conclusively demonstrate the existence of a torebody 
drag— base drag optimal point. The data demonstrate that 
the base drag coefficient corresponding to the drag 
minimum lies between 0.225 and 0.275, referenced to 
the base area. Most importantly, the data show a drag 
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reduction of approximately 15 percent when the drag 
optimum is reached. When this drag reduction is scaled 
to the X-33 base area, drag savings approaching 
45,000 N (10,000 lbf) can be realized. 

Nomenclature 

Acronyms 

CFD computational fluid dynamics 

LASRE Linear Aerospike SR-7 1 Experiment 

Symbols 

A slope parameter 

a^a ,, a a 3 forebody pressure distribution curve-fit 
coefficients 

law-of-the-wake bias parameter 
model span, cm 

base pressure distribution curve-fit 
coefficients 

intercept parameter 

drag coefficient 

base pressure drag coefficient 

forebody pressure drag coefficient 

zero-lift free-stream total drag coefficient 

“viscous” forebody drag coefficient 

local skin-friction coefficient 
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Phase 
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Ps 
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pressure coefficient 

average base pressure coefficient 

average forebody pressure coefficient 

section drag, N/m 

longitudinal pressure gradient on 
model, kPa/m 

longitudinal gradient of the boundary- 
layer momentum thickness 

expectation operator 

wake or boundary-layer shape parameter, 

H = 5/0 

base height, cm 

measurement index 

model length, cm 

number of repeated pressure scans 

number of data points 

local static pressure, kPa 

free-stream static pressure ratio in wind 
tunnel 

free-stream static pressure ahead of 
wind-tunnel model, kPa 


x 


y 



axial location within wind tunnel, cm 

/' lh scalar component of independent 
variable vector 

lateral coordinate (for wake, boundary 
layer, or base area), cm 

nondimensional boundary-layer 
coordinate 

output vector 

measurement vector 

/ th scalar component of measurement 
vector 

Clauser pressure gradient parameter 

friction velocity 

local curve-fit error for velocity 
distribution 

first variation of momentum 
thickness, cm 

gradient with respect to 5 

wake half-width, local boundary-layer 
thickness, cm 

initial estimate of wake half-width or 
boundary-layer thickness 


q dynamic pressure, kPa 

q ratu>( x } dynamic pressure ratio in wind tunnel 

q ^ free-stream dynamic pressure ahead of 

wind-tunnel model, kPa 

Re^ Reynolds number based on model 

length L 

Re v Reynolds number based on local axial 

coordinate x 

r leading-edge radius, cm 

U velocity at the edge of the wake or 

boundary layer, m/sec 

U ^ free-stream velocity ahead of the 

wind-tunnel model, m/sec 

u mi/i minimum velocity in wake velocity 

profile, m/sec 

«(y) local velocity distribution (in wake or 

boundary layer), m/sec 

u* nondimensional boundary-layer velocity 

X independent variable vector 


5* boundary-layer displacement 

thickness, cm 

5* u . wake displacement thickness, cm 

£ dummy integration variable 

9 forebody surface incidence angle, deg 

9 u . wake momentum thickness, cm 

9^ free-stream momentum thickness, cm 

K law-of-the-wake slope parameter 

K v equivalent sand-grain roughness of 

surface, cm 

K-e energy dissipation 

A Dirac delta function 

X roughness overlay ‘land” thickness, cm 

p sample mean 

E, dummy integration variable 

n wake parameter 

p air density, kg/cm 2 3 
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I 

roughness overlay slot thickness, cm 

2 

® base 

curve-fit squared error for base ports 

2 

® side 

curve-fit squared error for side ports 

2 

a e 

sample variance for pressure port 
incidence angle 0 

X 

roughness overlay shim thickness, cm 

i 

V'c 

L £>, 

base 

mean-square error in base drag 
coefficient estimate 

V 2 C 

^ forebody 

mean-square error in forebody pressure 
coefficient 

V 2 C 

C n 

u o 

mean-square error in viscous forebody 
drag coefficient estimate 

2 

V r 

mean-square error in viscous forebody 
drag coefficient estimate 

^AO 

mean-square error in momentum 
thickness estimate, crrr 

\p 2 

T A u/U e 

mean-square error in velocity profile 
curve fit, cm 2 

SuoerscriDts. Subscripts, and Mathematical Operators 

/ 

measurement index 

j 

pressure port index 

(*) 

iteration index 


estimated parameter 

A 

variational operator 

T 

vector transpose 


Introduction 


Designs advocated for the current generation of 
reusable launch or space-access vehicles are derived 
from variations of the original lifting-body concept. 1 
For many reasons, these designs all have large base 
areas compared with those of conventional aircraft. For 
example, the large base areas of the X-33 and Venture 
Star configurations are required to accommodate the 
aerospike rocket engines. The base area is highly 
separated, resulting in large negative base pressure 
coefficients. Because of the large base-to-wetted-area 
ratios of these vehicles, the base drag comprises the 
majority of the overall vehicle drag. The resulting low 
lift-to-drag ratios result in very steep approach glide 
slopes. These steep approach angles present difficult 
energy management tasks for autonomous reentry 
systems. Any decrease in base drag potentially can 
significantly improve the overall vehicle performance 


and make the autonomous reentry and landing task less 
difficult. 

An early body of experimental work conducted in the 
late 1950's and early 1960’s by Hoerner 2 offers a 
potential solution to the reusable launch vehicle (RLV) 
base drag problem. For blunt-based objects with heavily 
separated base areas, a correlation between the base 
pressure drag and the “viscous” forebody drag has been 
demonstrated. This paper presents the results of a series 
of wind-tunnel experiments that exploit this forebody- 
to-base drag relationship to reduce the overall drag of a 
simple blunt-based configuration by adding precise 
levels of roughness to the forebody. 

Background 

For blunt-based objects whose base areas are heavily 
separated, a clear relationship between base drag and 
the viscous forebody drag has been demonstrated by 
Hoerner. 2 In this paper, the viscous forebody drag is 
defined as the axial projection of the integral of all 
viscous forces acting on the vehicle forebody. These 
viscous forces include surface skin friction, frictional 
effects of forebody flow separation, and parasite drag. 
Axial forces resulting from the forebody pressure 
distribution are considered separately from the viscous 
forebody drag in this paper. 

Figure 1 shows subsonic drag data taken from 
Hoerner for two- and three-dimensional projectiles. The 
three-dimensional curve fit ot the data was originally 
published by Hoerner. 2 The two-dimensional curve fit is 
a new fit of Hoerner’s original data. The authors of this 
paper believe that this new fit is a better representation 
of the base drag data. 

An important feature is the trend for decreasing base 
drag as the viscous forebody drag increases (fig. 1). This 
base drag reduction is a result of boundary-layer effects 
at the vehicle base. The surface boundary layer acts as 
an insulator between the external flow and the separated 
air behind the base. As the forebody drag increases, the 
boundary-layer thickness at the forebody aft also 
increases. This increase reduces the effectiveness of the 
“jet pump” caused by the shearing ot the external flow 
on the separated flow behind the base region. 

Vehicle configurations w'ith large base drag 
coefficients lie on the steep portion of Hoerner’s curve, 
where a small increment in the forebody friction drag 
should result in a relatively large decrease in the base 
drag. Conceptually, if the added increment in viscous 
forebody drag is optimized with respect to the base 
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drag, then reducing the overall drag of the configuration 
may be possible. Figure 2 shows this drag optimization, 
based on curve fits of Hoemer's data. These data clearly 
illustrate the concept of the “drag bucket;’ 

Another important feature of the data shown in figures 
1 and 2 is that for the same viscous forebody drag, two- 
dimensional objects tend to have a significantly larger 
base drag than three-dimensional objects. These drag 
differences result from periodic shedding in the base 
region where a von Karman vortex street structure 3, 4 of 
evenly spaced vortices of alternating strengths sets up 
within the wake. In general, the base flow around 
three-dimensional objects is characterized by very- 
broadband (frequency) flow disturbances; the periodic 
flow phenomenon is far less pronounced than for two- 
dimensional objects. The base pressure under 

nonperiodic (three-dimensional) flow conditions is 
considerably higher (equating to lower base drag) than 
under similar conditions in a periodic (two-dimensional) 
flow. 

The ramifications of this two-dimensional-three- 
dimensional base drag difference become 

extremely important when one considers full-scale, 
high-Reynolds number flight vehicles. Saltzman, et al. 4 5 
have compiled subsonic drag data from vehicles 
configured for hypersonic flight. This compendium 
includes flight data for the X-15, M2-F1, M2-F3, 
X-24A, and X-24B vehicles; the Space Shuttle; and the 
Linear Aerospike SR-71 Experiment (LASRE). These 
data are compared to the two- and three-dimensional 
mathematical models derived from Hoerner’s data 
(fig. 3). The full-scale flight data clearly agree more 
closely with the two-dimensional curve than the three- 
dimensional one. For full-scale configurations, the flow 
appears to be locally two-dimensional and allows the 
trailing vortex street to become well-established. 
Figure 4 shows direct visual proof of this assertion, as a 
periodic vortex structure is clearly visible trailing 
behind the M2-F1 vehicle. 

The data shown in figures 1-3 imply that large-scale, 
blunt-based vehicles are quasi-two-dimensional, and 
configurations with a base drag coefficient greater than 
approximately 0.30 (referenced to the base area) will lie 
on the left side of Hoerner’s curve. These configurations 
may be considered to be suboptimal with respect to the 
viscous forebody drag coefficient. Incrementally 
increasing the viscous forebody drag theoretically 
should lower the overall drag of the configuration. 


The Linear Aerospike SR-71 Experiment 

Flight test results from the LASRE drag reduction 
experiment 6 provide some incomplete validation of the 
above hypothesis. The LASRE was a flight test of an 
approximately 20-percent half-span model of an X-33 
forebody model mounted on top of the NASA SR-71 
aircraft. The LASRE sought to reduce base drag by 
adding a small amount of surface roughness to the 
model forebody. The model was instrumented with load 
cells that allowed a six-degree-of-freedom measurement 
of forces and moments, and with surface pressure ports 
that allowed the model forebody pressure and base drag 
to be numerically integrated. 

The LASRE verified that the base drag was reduced 
by as much as 15 percent; unfortunately, the overall drag 
of the configuration was not reduced. The methods for 
applying the forebody sand-grain roughness were 
believed to be too crude to achieve an overall drag 
reduction. Further tests under a more controlled flow’ 
environment were clearly required. 

Wind-Tunnel Tests 

A series of low-speed, two-dimensional, wind-tunnel 
tests was conducted to study the potential for 
minimizing the total configuration drag using surface 
roughness increments. In these tests, a leading-edge 
cylinder with a blunt afterbody was tested. The full- 
scale flight data (figs. 3—4) demonstrate that the results 
of the two-dimensional tests should be generally 
applicable to large-scale, three-dimensional vehicles. In 
fact, with regard to the comparisons shown in figure 3, 
tests performed using two-dimensional models were 
believed to be more representative of the large-scale 
flight vehicles than those performed with three- 
dimensional models. The series of tests had two primary 
objectives: 

1. Test the hypothesis regarding forebody roughness 

in a systematic manner to conclusively 

demonstrate existence of a viscous forebody drag- 
base drag optimum (the “drag bucket”). 

2. Establish a criterion for when forebody drag is 
suboptimal (that is, at what point does increasing 
forebody drag result in an overall drag reduction). 

Wind-Tunnel Model Description 

Figure 5 shows a three-view drawing of the 
wind-tunnel model. The machined-aluminum model 
consists of a 2.54-cm-diameter (1-in.) cylindrical 
leading edge with a flat-sided afterbody 11.43-cm 
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(4.5-in.) long. Removable aluminum plates on the sides 
of the model allow various levels of surface roughness 
to be tested by interchanging the plates. The base-to- 
wetted area of the model is approximately 10.7 percent. 
Figure 6 shows the model mounted in the wind tunnel. 

The forebody roughness of the model was increased 
by bonding micromachined brass overlays to the side 
plates. Figure 5 shows a sample of this roughness 
“screen” overlaid on the top view of the model. These 
“screens” consist of a series of transverse bars with the 
shim (x), slot (I), and “land” (X) dimensions 
determining the roughness of the surface. Figure 7 
shows the geometric layout for these bar grid overlays. 
A single overlay geometry using lands and slots aligned 
parallel to the direction of flow was also tested. Table 1 
shows the geometries tested, and the equivalent surface 
roughness (K v ) derived from empirical-fit formulae 
presented in Mills. 7 


Table 1. Screen overlay roughness dimensions. 


Configuration 

number 

X , cm 

Z . cm 

x , cm 

K y ,cm 

l 

0.0000 

0.0000 

0.0000 

0.0000* 

2 

0.005 1 

0.0051 

0 0051 

0.0163** 

3 

0.0254 

0.0381 

00254 

0.1143 

4 

0.0508 

0.1016 

0.0508 

0.2896 

5 

0.0508 

0.2032 

0.0508 

0.4854 

6 

0.1016 

0.2540 

0 1016 

0.6911 


* Smooth model 
** Parallel bars 


Wind-Tunnel Description 

The model was tested in a low-speed wind tunnel at 
the NASA Dryden Flight Research Center (Edwards, 
California). The ambient, open-cycle tunnel has a test 
section approximately 10 by 25 cm (4 by 10 in.). An 
alternating current (A/C) motor uses a squirrel-cage fan 
located at the downstream end to pull air through the 
tunnel. When the model was mounted in the tunnel test 
section, the total blockage was 10 percent. This level of 
blockage is considered high for traditional wind-tunnel 
testing. 

The primary effect of the blockage was to accelerate 
the flow around the model forebody, causing a rise in the 
dynamic pressure and a drop in the static pressure along 
the sides of the tunnel wall (outside of the tunnel wall 
boundary layer). The dynamic pressure rise (static 
pressure drop) was taken into account by calibrating 


local total and static pressure ratios — referenced to the 
dynamic and static pressure ahead of the model — as a 
function of the axial position in the tunnel. Figure 8 
shows this calibration plot. At each pressure 
measurement location, the derived dynamic pressure 
was used to compute the local pressure coefficient. 


C p (x) = 


[P^ ~ Psra,J X)P ~_ 


QratioWq* 


( 1 ) 


With the model mounted in the wind tunnel, a 
maximum free-stream airspeed of approximately 
28.0 m/sec (92 ft/sec) was achieved. Based on the model 
length, this free-stream velocity translates to a Reynolds 
number (Re^) of approximately 2.25 x 10 .Tests were 
also performed at airspeeds of approximately 
14.6 m/sec (48 ft/sec). The corresponding Re, for these 
lower-speed tests was approximately 1.25 x 10 . The 
wind-tunnel turbulence intensity levels were sufficiently 
large that the model flow was turbulent beginning at the 
leading edge. 


Instrumentation 

All test measurements were performed using only 
pressure instrumentation. The methods used to interpret 
the measurements are presented in the "Analysis 
Methods” section. The tunnel itself was instrumented 
with series of static pressure taps along the side of the 
tunnel. Total (reference) pressure levels were sensed 
with a pitot probe placed five model lengths ahead of the 
model. A total of 16 pressure taps was distributed 
around the centerline of the model: 5 ports on the model 
forebody, 8 ports placed along the sides of the model, 
and 3 ports placed on the model base. These port 
locations allowed body pressure forces to be accurately 
integrated. Figure 5 shows the locations of the 16 model 
pressure ports. Several leading-edge ports can be seen 
on the model mounted in the tunnel (fig. 6). 

The total model drag coefficient was measured by 
wake velocity profiles sensed using a traversing pitot- 
static probe. Both local total and static pressures were 
sensed by this probe. The probe tip was placed 12.7 cm 
(5 in.) aft of the model base area. The wake probe tip 
diameter was approximately 0.025 cm. Similar 
momentum-defect measurements for skin friction were 
performed at the model aft using a traversing 
boundary- layer pitot probe. For the boundary-layer 
profiles, only local total pressure was measured by the 
traversing probe. Local static pressure was assumed 
constant across the depth of the boundary layer. 
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Free-stream static pressure at the model base was sensed 
by a side port on the tunnel wall. The boundary-layer 
probe tip diameter was approximately 0.02 cm. Figure 6 
shows the wake and boundary-layer probes mounted in 
the tunnel. The probe positions relative to the centerline 
of the model were measured using a digital micrometer. 
The estimated accuracy of the digital positioning sensor 
was approximately 0.0025 cm (0.001 in.). 

All of the model, tunnel wall, and traversing probe 
pressure data were sensed with a highly accurate set of 
digital (RS-422) scanning pressure modules. These data 
were recorded by a laptop computer using the serial port 
to perform individual channel addressing. Full-scale 
span of these differential pressure modules was 
±2.490 kPa (±52.0 lbf/fr). The manufacturer’s accuracy 
specifications for the differential pressure measurements 
is ±0.05 percent of full scale, or approximately 
±0.00125 kPa (±0.026 Ibf/ft 2 ). The differential pressure 
transducers were referenced to the pitot probe placed 
approximately 64 cm (25 in.) ahead of the model. The 
reference pitot pressure was sensed with a highly 
accurate absolute pressure manometer. The estimated 
accuracy for the absolute reference pressure 
measurement is approximately ±0.010 kPa 
(±0. 16 lbf/fr). The reference temperature was sensed 
externally to the tunnel using a type “T” thermocouple 
with an estimated accuracy of approximately 
±0.5 °C (±0.9 °F). 

Test Procedures 

The low dynamic pressure levels — less than 
0.4788 kPa (10 lbf/fr) — during this series of wind- 
tunnel tests required that data be taken with great 
consistency to minimize the effects of experimental 
procedure on the overall errors. For all test conditions 
and configurations, the transducers were zeroed prior to 
testing, and the model angle of attack was set to zero by 
comparison of the left and right surface model 
pressures. To set the zero angle-of-attack position, the 
model position was perturbed until the left and right 
surface pressure curves lay directly on top of each other. 

Transducer Zeroing 

Although the electronically scanned pressure 
transducers have a built-in feature that allows the 
transducers to be zeroed on-line, experimentation 
determined that a superior level of bias correction was 
achieved when the transducers were manually zeroed 
before each data run. Transducer biases were evaluated 
by taking readings with the tunnel in the “off” position 
(zero airspeed). In this zeroing process, each pressure 


port was addressed a total of 100 times and these data 
samples were averaged to minimize the effects of 
random sensor errors. The resulting zero readings were 
written to an archival file for later use by postprocessing 
analysis algorithms. 

Surface Pressure Scans 

The pressure scans read data from the 16 model 
pressure ports as well as the total and static pressure 
levels in the tunnel. For each configuration tested — that 
is, each different grid pattern or airspeed — the pressure 
scans were repeated ten times. For each of the ten 
measurement sequences, the zeroing procedure was 
performed and the tunnel was activated and allowed to 
stabilize. Typically, 100 individual data samples were 
averaged for each data run to minimize the effects of 
random measurement errors and tunnel turbulence. 
After ten pressure scans were taken for each 
configuration, the data were converted to pressure 
coefficients by postprocessing algorithms and the 
pressure coefficients data were averaged. The standard 
deviation of the ten measurement sequences data was 
used as a representation of the end-to-end accuracy of 
the measurement system. Typically the end-to-end 
pressure coefficient error varied between ±0.003 
and ±0.005. 

Wake and Boundary-Layer Surveys 

For the wake surveys, each data point consists of a 
pitot and a static-pressure measurement taken at a single 
lateral offset (y) from the model centerline. For the 
boundary-layer surveys, each data point consists of a 
pitot measurement taken at a lateral offset and a wall 
static pressure measurement. For each data point, 100 
data samples were averaged to minimize the effects of 
random measurement errors and tunnel turbulence. To 
completely define the wake profile, approximately 200 
y-position data points were required. For early tests in 
the tunnel, the entire wake profile was measured. These 
data were so symmetrically distributed that as a time- 
saving measure, later tests only surveyed one-half of the 
wake profile. 

Because of the large number of data samples 
(approximately 20,000) required to define the wake for 
each measurement configuration, completing each of 
the wake surveys ten times as was done with the 
pressure survey data was considered impractical. 
Instead, each wake survey was performed twice and the 
resulting data were interleaved to form a single local 
velocity distribution profile. At the beginning of each of 
the two wake surveys, the probe sensor zero readings 
were taken and written to an archival file for use by the 
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postprocessing routines. When computed, transducer 
biases were assumed constant for the duration of each 
wake survey. 


suboptimal side of Hoemer's drag curve. Thus, by 
adding roughness to the forebody, the overall drag 
coefficient should be reduced. 


Analysis Methods 


Wake Profile Analysis 


This section derives the analysis methods used in this 
series of wind-tunnel tests. A baseline set of two- 
dimensional, incompressible, computational fluid 
dynamics (CFD) calculations will be presented first. 
Next, the viscous calculations used to convert the 
measured wind-tunnel pressures data into the various 
components of the drag coefficient will be presented. 
For each analysis method presented in this section, an 
error analysis is also presented in the appendix. 

Computational Fluid Dyna mics Analysis 

The CFD calculations were performed to give pretest 
drag predictions to verify that the smooth model 
configuration lay on the suboptimal portion of Hoerner s 
base drag curve. 

Only the CFD estimates of forebody pressure and 
base drag coefficients were used for the pretest drag 
predictions. Not enough computational cells were 
embedded within the boundary layer to allow the 
skin-friction coefficient to be accurately computed, 
using the CFD data. The integrated skin drag coefficient 
was predicted using the two-dimensional Floerner drag 
model. 


This analysis method fits the wind-tunnel wake data 
with a symmetric ik cosine law” velocity distribution 
profile of the form 


»(y) _ I 

u m 2 




COSI Tt= 

^ o 


( 2 ) 


In equation (2), u min is the minimum velocity in the 
wake, y is the lateral distance outwards from the center 
of the wake, U is the velocity at the edge of the wake, 
u(y) is the local velocity within the wake, and 5 is the 
wake half-width. A least-squares method was used to 
curve-fit the measured velocity distribution data to the 
profile assumed in equation (2). In this method, equation 
(2) is rewritten as a linear system of the form 

Z (mrav) = AX (k) + C (3) 


where 


The CFD flow calculations were performed using a 
commercially available code.** The core solver for this 
code features a finite-volume, cell-centered 
discretization, and uses a time-accurate, ”PISO” 
(pressure-implicit with splitting of operators) solution 
algorithm to solve the integral form of the 
Navier-Stokes equations. Although the code has 
compressibility and transient solution capabilities, only 
the incompressible steady-state solution was used in this 
analysis. The analysis was set up to force turbulence at 
the leading edge of the model. For this analysis, a 
simple K- £ (energy-dissipation) turbulence model was 
used. 

Figure 9 shows the predicted CFD model flow field. 
The CFD solutions clearly show a periodic vortex 
structure trailing the model. When the pressure forces 
are summed along the surface of the model and 
projected perpendicular to the longitudinal axis, the 
integrated forebody pressure coefficient is 
approximately —0.018 and the integrated base drag 
coefficient is approximately 0.035. Based on data shown 
in figure 2, the smooth model should lie on the 
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A simple least-squirms method is used to solve for 
estimates of the slope and intercept parameters, A 
and C: 
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(A) 


v m 


= A lk) + t w 


(4) 


Using a first-order perturbation, equation (4) can be 
“updated’' using nonlinear regression to get a refined 
value for 5 : 


7 (meas) ^(*) 


<*>* 


-Z = A VcX 


,(k) 


S a + ,) -8 (i) 


(5) 


and the cosine velocity distribution law gives a 
reasonable curve fit. Note that the center of the wake 
appears to contain a significant amount of turbulence 
that significantly decreases near the edge of the wake. 

When the velocity profile has been curve-fit, 
equation (1) is substituted into the equations for the 
wake displacement and momentum thickness and 
analytically evaluated to give 


where 


( \ 



After extensive algebra, the least-squares solution to 
equations (5) and (6) can be written as 
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In equations (8) and (9), u(y) is the local velocity in 
the wake at lateral offset location y , and U is the local 
velocity at the edge of the wake. The free-stream 
momentum thickness is calculated from the local 
momentum thickness using the well-known Squire- 
Young formula: 10 


5 (* + U S (A) 

o =0 
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Assuming that a starting value for the wake 
half-width, is known beforehand (from visual 

inspection of the wake data), equations (4)— (7) are 
solved iteratively until convergence. Convergence 
typically takes less than ten iterations. 

Figure 10 shows an example wake curve fit compared 
with the wind-tunnel data. These data were obtained 
from the smooth model configuration tested at 
Rej = 2.25x10 . The turbulent wake extends beyond 
the lateral boundaries of the wind-tunnel model by 
approximately 3 cm. The wake structure is symmetric 


= 0. 


i/ 


i[H + 5| 


(10) 


Equation (10) corrects for the effects of the wind 
tunnel blockage described earlier in this paper. In 
equation (10), H is the wake shape parameter defined by 

8 * h . 

H = ^ (11) 


The free-stream drag coefficient is computed from the 
normalized section drag 


C n - 


D' 




l ba\e 


1 base 


( 12 ) 


An approximate accounting of overall error in the 
wake drag coefficient can be performed using a linear 
perturbation analysis. The appendix shows this 
linearized error analysis. 
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Boundary-Layer Profile Analysis 

The forebody skin friction coefficient is evaluated 
using the boundary-layer velocity profiles in a similar 
manner as the wake analysis presented earlier. In this 
case, however. Coles’ “law of the wake,” 

ln[ y + ] + 2U sin" 



ny 

25. 


+ B 


(13) 


is curve-fit to the local velocity profile data. The law of 
the wake is a very general experimental correlation for 
turbulent boundary layers, and relates the 
nondimensional velocity 


«oo 



(14) 


to the nondimensionalized boundary-layer coordinate 


+ 
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(15) 


dP /dx is the longitudinal pressure gradient at the 
edge of the boundary layer. Based on the correlation of 
equation (17), the numerical value of II corresponding 
to a zero pressure gradient flow is approximately 0.426. 
Earlier authors have placed this zero gradient value at 
approximately 0.5 10 and 0.55 7 For this analysis, the 
more modern value recommended by Das is used. A 
value for n greater than the zero gradient value (0.426) 
corresponds to an adverse pressure gradient. A value for 
n less than the zero gradient value (0.426) corresponds 
to a favorable pressure gradient. 10 

Following the procedure used earlier with the wake 
integral analysis, equation (16) is rewritten as a linear 
system of the form 


C 1 


v i 


= r 




x 



_ n _ 


In equations (13)-(15), 5 is local boundary-layer 
thickness, k is the iaw-of-the-wake slope parameter, B 
is the law-of-the-wake bias parameter, n is the wake 
pressure gradient parameter, and cy is the local 
skin-friction coefficient. The accepted “best value for 
k currently is 0.41. 10 The bias parameter, varies with 
the level of surface roughness and for a smooth plate has 
a numerical value of approximately 5.0. Re v is the 
Reynolds number based on the local axial coordinate, x. 

The roughness dependent bias term can be eliminated 
from equation (13) by expressing the law of the wake in 
terms of the local “velocity detect : 


where 



In equations (18M19), the subscript i is the 
measurement, and the superscript (k) is the iteration 
index. After some extensive algebra, the least-squares 
solution to equation (19) can be written as 



2 TT cos 




(16) 


The wake parameter, FI , is proportional to the local 
longitudinal pressure gradient. Das 10- has established 
an empirical correlation that relates the wake parameter 
to the more familiar “Clauser parameter,” 17 p , where 
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In equation (17), 5* is the local displacement 
thickness, c . is the local skin-friction coefficient, and 
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( 20 ) 


Using a first-order perturbation with respect to 5, 
equation (20) can be updated using nonlinear regression 
to get a refined value for 5 : 
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The resulting updated equation for 5 is 
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When the variational algorithm of equations 
(21 )— ( 23 ) is modified to allow direct estimation of IT 
along with 8 and Cy , the equations rapidly diverge. To 
circumvent this numerical problem, H was selected for 
this analysis to give the best overall fit consistency. This 
procedure typically consisted of selecting a starting 

value for n and then computing c f and 5 by 

’ * 

iteratively solving equations (21 )— ( 23) until 

convergence. At this point, the value for n was varied 
by a small amount and the iterative algorithm was 
repeated. If the total fit error improved, then FI was 
again varied in the same direction; if not, then the value 
was varied in the opposite direction. Using this ad hoc 
procedure, a minimum tit error is typically reached after 
less than ten trials. 


n = 1.032; and a l/7 lh -power-curve exponential curve 
fit. Analysis of equation (17) presented in White 10 
shows that n= 1.032 corresponds to a weak adverse 
pressure gradient. The model data presented in the 
‘"Results and Discussion” section support this 
conclusion. Clearly, the curve fit using n= 1.032 gives 
overall fit consistency. 

The estimated values for 5 , cy , and n are used to 
calculate the local momentum and displacement 
thickness by integrating the law of the wake across the 
depth of the boundary layer. As derived in White, 10 the 
resulting expressions for the displacement and 
momentum thickness are 


5 * _ Ff, i + n 

6 - V 2 k 


and 


(24) 


e 

8 



2 + 3.2FI+ 1 .5 n 2 ] 


(25) 


For simplicity, the effect of the local laminar sublayer 
is ignored in equations (24) and (25). For the Reynolds 
numbers tested, earlier analysis estimates that ignoring 
the laminar sublayer introduces integral errors of less 
than 0.2 percent. When the local momentum and 
displacement thickness have been evaluated, then the 
integrated viscous forebody drag coefficient can be 
evaluated using the “Clauser” form of the von Karman 
momentum equation, 10 


dQ 0 + H^ Cf ' - Cf ' 

Zr {2 + H) H T - T 


(26) 


In equation (26), H = 5*/0 is the boundary-layer 
shape parameter. The Clauser parameter, (3 , is related to 
the local pressure gradient, the displacement thickness, 
and the local skin-friction coefficient as 


Figure 1 1 shows an example boundary-layer curve fit 
compared with the wind-tunnel data. The normalized 
velocity distribution is plotted against the normalized 
position within the boundary layer. These data were 

obtained from the smooth model configuration tested at 

5 

Re^ = 2.25 x 10 . Three fit curves are plotted here: a 
law-of-the-wake curve fit with n = 0.426 (zero pressure 
gradient); a law-of-the-wake curve fit with the wake 
parameter adjusted to give the minimum fit error. 


p 





5* dP e 



2 dx 


(27) 


Solving equations (26) and (27) for the local skin- 
friction coefficient gives 


= H 

C f< ~dx\ll + (2 + W)P! 


( 28 ) 
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As demonstrated by Clauser, 13 for small-to-moderate 
pressure gradients, the terms on the right side of 

equation (28), — — — — tttttt , are approximately 

[H + (2 + n )pj 

constant. Integrating equation (28) along the forebody 
length, L , gives 


for the taper of the base pressure near the outer edges of 
the model. In this curve-fitting scheme, ports 7 and 1 1 
were weighted one-half as much as the three base area 
ports (ports 8, 9, and 10). This weighting scheme was 
selected to give a base drag taper correction factor of 
approximately 0.925. This correction factor is suggested 
by Saltzman, et al. 5 for full-scale flight vehicles. 


C F 


- ifVjg H 

dx iN + {2 + H)$] 
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(29) 


As with the earlier wake analysis, an approximate 
accounting of overall error in the wake drag coefficient 
can be performed using a linear perturbation analysis. 
The appendix shows this linearized error analysis. 

Forebodv Pressure Analysis 

The forebody pressure coefficient was evaluated by 
curve-fitting the pressure distributions as a function of 
local incidence angle. 0 . For the forebody data, seven 
forebody pressures — ports 1, 2, 3, 4, 14, 15, and 16 
(fig. 5) — are curve-fit with a third-order polynomial. 
The forebodv pressure drag coefficient is analytically 
given by the surface integral 
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= « () + 0.5708 a , + 0.4674 a 2 +0.4510 


Figure 12 shows a plot of a sample forebody curve fit. 
These smooth model data were measured with the wind 
tunnel operating at Re^ = 2.25 x 10 . The forebody 
pressure coefficient data is plotted as a function of the 
local incidence angle. The upper and lower surface 
pressure data lie nearly superimposed on each other, so 
not surprisingly, the third-order curve-fit closely 
matches the pressure coefficient data. 

Base Pressure Analysis 

The base pressure coefficient was evaluated by 
curve-fitting the base pressure distributions as a function 
of the lateral offset coordinate, y . For the base pressure 
data, five base area pressure ports — ports 7, 8, 9, 10, and 
11 (fig. 5) — were curve-fit with a fourth-order 
polynomial. The pressure ports on the sides of the model 
(ports 7 and 1 1) were included in the curve fit to account 


The base pressure drag coefficient is given 
analytically by the evaluating the surface integral 





C p [y] dy = 



dy 


( 31 ) 


= b 0 + 0.0833 b 2 + 0.0125 b 4 


Figure 13 shows a sample base pressure distribution 
curve fit. These data were measured on the smooth 
model with the wind tunnel operating at an approximate 
Reynolds number of 2.25 x 10 , based on model length. 

Results and Discussion 

The wind-tunnel data clearly support the earlier CFD 
predictions that the smooth model will lie on the 
suboptimal side of Hoerner’s curve. The suboptimal 
hypothesis is most clearly demonstrated by examining 
the base area pressure distributions. Figure 14 shows 
these results. The base pressure coefficients are plotted 
here as a function of v for various surface grid patterns. 
Figure 14(a) shows the pressure distributions for 
Re L = 2.25 x 10 5 , and figure 14(b) shows the pressure 
distributions for R e ^ = 1.25 xlO 5 - interestingly, the 
surface pattern with fine-mesh parallel slots and lands 
causes the base drag to dramatically rise (and have 
lower base pressure coefficients) when compared with 
the smooth surface model. Conversely, the surface 
pattern with transverse slots and lands causes the base 
drag to gradually lower (and have higher base pressure 
coefficients) when compared to the smooth surface 
model. 

A similar behavior was observed by Krishnan, et al., 14 
when the authors added riblet 1 ^ structures to the 
forebody of an axisymmetric w ind-tunnel model w ith a 
blunt base. The authors’ intents were that the riblets 
would lower base drag; however, the results were 
opposite of expectations. When Krishnan’s results and 
the data presented in figure 14 are interpreted 
considering Hoerner’s curve (fig. 3), the rising base drag 
is completely reasonable. The grid pattern with parallel 
slots and lands has the effect of acting like riblets on the 
model forebody. The riblet structures have the effect of 
lowering the forebody drag coefficient. Because the 
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forebody skin drag coefficient is lowered, the base drag 
is expected to correspondingly increase. Clearly, riblets 
should not be used in conjunction with “suboptimal” 
configurations that have highly separated base regions; 
their effect will cause the base drag to rise. 

Figure 15 shows results from the wind-tunnel tests 
that further illustrate this concept. The measured base 
drag coefficient is plotted with the viscous forebody 
drag coefficient calculated from the boundary-layer 
survey data. These data are compared to the curve fit of 
Hoerner's two-dimensional data from figure 1 -The open 
symbols represent data for Re^ = 2.25 x 10 and th^e 
closed symbols represent data for Re^ = 1.25 x 10 . 
The error bars show the expected “l-O* standard 
deviations based on the error analyses presented earlier. 
The agreement with the curve fit of Hoerner s data is 
reasonably good. 

Figure 16 shows the model total drag coefficient data 
(as calculated from the wake survey data) plotted with 
the viscous forebody drag coefficient. The error bars 
show the expected “1-a” standard deviations based on 
the error analyses presented earlier. Figure 16 also 
shows the predicted drag curve defined using Hoerner’s 
two-dimensional curve from figure 1, the viscous 
forebody drag measurement, and the model forebody 
drag coefficient predicted (-0.018) by the CFD 
solutions. Note that, with the exception of the data for 
the parallel grid (riblets) overlay, the agreement with the 
predicted drag curve is very good. 

The disagreement for the parallel grid test points is 
caused by a sharp rise in the forebody pressure 
coefficients. Figure 17 shows these data. The forebody 
pressure distributions for all of the grids are plotted here 
as a function of the local incidence angle. Figure 17(a) 
shows the higher Reynolds number data (2.25 x 10 ), 
and figure 17(b) shows lower Reynolds number data 
( 1.25 x 10^). The transverse grid patterns do not 
significantly alter the forebody pressure distribution; 
however, the forebody pressure data are considerably 
higher for the parallel grid pattern. The parallel grid data 
are clearly an anomaly. The reasons for this pressure 
anomaly are not clear at this point, but the parallel grid 
possibly caused relaminarization of the flow and 
induced a localized separation. This anomaly requires 
further investigation. 

Most importantly, the data shown in figure 16 
demonstrate the existence of a drag minimum with 
regard to the viscous forebody drag coefficient. The 


elusive “drag bucket" is clearly defined and the primary 
hypothesis of this paper is conclusively proven. The 
drag reduction from the smooth model configuration to 
the optimum point is approximately 15 percent. Also, 
comparison of figure 15 with figure 16 shows that the 
base drag coefficient corresponding to the total drag 
coefficient minimum lies somewhere between 0.225 and 
0.275. This value is a bit lower than the 0.25-0.30 range 
predicted by analysis of Hoemer’s original data (figs. 2 
and 3). 

Summary and Concluding Remarks 

Current designs of transatmospheric crew return and 
reusable launch vehicles have extremely large base-to- 
wetted area ratios when compared to conventional 
vehicle designs. These truncated base areas are highly 
separated, resulting in large, negative, base pressure 
coefficients. Because of the large base-to-wetted-area 
ratio, base drag makes up the majority of overall vehicle 
drag. Any reduction in base drag directly improves 
vehicle performance, resulting in an enhanced lift-to- 
drag ratio, extended range, and a less-severe approach 
glide slope. 

Early work performed on blunt-based bodies offers a 
potential solution. For blunt-based bodies, a direct 
correlation exists between base and “viscous” forebody 
drag. As the forebody drag coefficient increases, the 
base drag of the projectile generally tends to decrease. 
This base drag reduction results from boundary-layer 
effects at the vehicle base. Conceptually, if the added 
increment in forebody skin drag is optimized with 
respect to the base drag reduction, then reducing the 
overall drag of the configuration may be possible. 

In order to test the above concept, a series of 
small-scale wind-tunnel tests was conducted. In these 
tests, a two-dimensional cylinder with a blunt afterbody 
was tested. The series of tests had two primary 
objectives: to test the forebody roughness hypothesis in 
a systematic manner to conclusively demonstrate 
existence of a “drag bucket”; and to establish a criterion 
for when forebody drag is suboptimal (that is, when will 
increasing forebody drag result in an overall drag 
reduction). 

This paper presents the wind-tunnel test results. Both 
primary objectives were satisfied. These wind-tunnel 
results conclusively demonstrate existence of a 
forebody drag optimum. Also, the wind-tunnel data 
demonstrate that the base drag coefficient corresponding 
to the total drag minimum lies somewhere between 
0.225 and 0.275. This optimality point is slightly lower 
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than the 0.25-0.30 range predicted by analysis of 
Hoerners original data. The use of parallel grid lines 
that emulate the effects of riblet structures on bodies 
with highly separated base regions will likely cause the 
total drag of the configuration to rise. Most importantly, 
the data show a peak drag reduction was approximately 
15 percent. When this 15-percent drag reduction is 
scaled to the size of the X-33 vehicle, the drag savings 
approaches approximately 45,000 N (10,000 lbf). 

Clearly, this experiment should be repeated for 
different ranges of Reynolds number and aspect ratios to 
determine if the lower optimality point indicated by the 
data is real. The methods should also be demonstrated as 
being effective in the presence of induced drag. Practical 


implementation methods that allow for on-line adaptive 
modification of the forebody drag coefficient to seek the 
optimal point should be explored and developed. The 
limits of practical applicability for this technology are 
unknown at this point. This drag reduction technology is 
still in its infancy; however, a wide spectra of potential 
users exist, including the aerospace, automotive, ground 
transport, and shipping industries. Use of this drag 
reduction technique offers the potential for decreased 
operating costs resulting from decreased overall fuel 
consumption. 


o Two-dimensional data (from Hoerner) 

Two-dimensional curve fit 

□ Three-dimensional data (from Hoerner) 
Three-dimensional curve fit 



Figure 1. The effect of the viscous forebody drag on the base drag of a blunt-based projectile. 
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Two-dimensional drag optimization 

Three-dimensional drag optimization 



Drag coefficients 

Base (two-dimensional model) 

Total (two-dimensional model) 

Base (three-dimensional model) 

Total (three-dimensional model) 

□ Base (flight data) 



000638 

Figure 3. Comparison of flight data to two- and three-dimensional drag models. 
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Figure 7. Bar grid surface overlay screen pattern. 
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1.3 


1.0002 
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Figure 9. Two-dimensional, incompressible CFD solutions of wind-tunnel model flow field (smooth configuration). 
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Cosine curve fit 

□ Measured data 



Figure 10. Example wind-tunnel wake data. 


Curve fit, n = 1.032 

Curve fit, FI = 0.426 

1/7th power law fit 

o Wind-tunnel test data 



Figure 1 1. Example wind-tunnel boundary-layer curve fit. 
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Figure 12. Example forebody pressure coefficient curve fit. 


Fourth-order 

polynomial curve fit 
O Wind-tunnel data points 



y* Cm 000649 


Figure 13. Example base pressure coefficient curve fit. 
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O Smooth model 

□ Transverse grid 2, K s - 0.1143 cm 
O Transverse grid 4, K g - 0.4854 cm 
A Parallel grid 1, K s - 0.0163 cm 
fcs. Transverse grid 3, K s ~ 0.2896 cm 



O Smooth model 

□ Transverse grid 2, K s - 0.1143 cm 
O Transverse grid 4, K s ~ 0.4854 cm 
A Parallel grid 1, K g - 0.0163 cm 
bs. Transverse grid 3, K s ~ 0.2896 cm 
Qi Transverse grid 5, K s - 0.6911 cm 



(b) Re, = 1.25 x 10 5 . 

Figure 14. Base pressure distributions for various grid patterns. 
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Figure 15. Comparison of wind-tunnel base drag data as a function of Hoerner’s curve. 
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Figure 16. Comparison of wind-tunnel wake survey data as a function of the predicted "drag bucket.” 


• Re L ~ 1.225 x10 s 
□ Re|_ - 2.250 x 10 5 



• Re L ~ 1.225x10 5 
□ Re L ~ 2.250 x 10 5 
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O Smooth model 

□ Transverse grid 2, K s - 0.1143 cm 
O Transverse grid 4, K s - 0.4854 cm 
A Parallel grid 1, K s - 0.0163 cm 
Transverse grid 3, K s ~ 0.2696 cm 
Di Transverse grid 5, k s ~ 0.6911 cm 



(a) Re L = 2.25 x 10 5 . 


O Smooth model 

□ Transverse grid 2, K s - 0.1143 cm 
O Transverse grid 4, K s - 0.4854 cm 
A Parallel grid 1, K s ~ 0.0163 cm 
Transverse grid 3, K s - 0.2896 cm 
IA Transverse grid 5, K s ~ 0.6911 cm 



(b) Re^ = 1.25 x 10 5 . 


Figure 17. Forebody pressure distributions for various grid patterns. 
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APPENDIX 

ERROR ANALYSIS METHODS 

Introduction 

An approximate accounting of overall error in the 
wind tunnel— derived drag coefficient estimates is 
derived herein. The wake error analysis is presented 
first; the boundary-layer skin-friction error analysis is 
presented next. The estimated errors in the forebody and 
base pressure drag coefficients are presented last. The 
error equations derived in this appendix were used to 
calculate the data point error bounds plotted in figures 
15 and 16. 

F.sti mating the Wake Drag Coefficient Errors 



An approximate accounting of overall error in the 
wake drag coefficient can be performed using a linear 
perturbation analysis. Using the fundamental definition 
for momentum thickness, linear perturbations can be 
expressed in terms of the velocity profile curve-fit error 

by taking the first variation of the forebody 
surface incidence angle, 0 , with respect to u(y)/ U e : 


A0 = A 


«(?)/!/, 






(A- 1 ) 


The mean-square error in momentum thickness is 
evaluated by taking the expectation of the square of 
equation (A-l). 


Assuming that random local errors in the velocity 
profile curve fit are uncorrelated, the expectation 
operation in equation (A-2) becomes 


E 
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(A-3) 


In equation (A-3), A^ is the Dirac delta function 16 

and v P 2 Am/( / is the mean-square error in the velocity 

distribution curve fit. Substituting equation (A-3) into 
equation (A-2), the interior integral reduces to the 
following: 



Substituting equation (A-4) into equation (A-2), and 
using the wake cosine law (eq. (2)) to evaluate the outer 
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integral, the approximate mean-square error in the 
momentum thickness reduces to the following: 
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A first-order perturbation of equation (A-5) gives the 
error equation for the free-stream momentum thickness 
estimate (assuming that errors in the velocity ratio and 
shape parameter are negligible when compared to the 
momentum thickness errors): 
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Finally, the approximate mean square error in total 
drag coefficient is 
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In equation (20), H is the shape parameter defined by 

equation (10), 5 is the boundary-layer thickness, and 
*> 

T‘a {4 /(f is the mean-square error in the velocity 
profile curve fit. 

Estimating the Forebody Viscous Drag 
Coefficient Errors 


The boundary-layer error analysis follows a nearly 
identical process when compared to the wake error 
analysis. The main exceptions are that the integrals are 
performed from {0,5} instead of {-5,5}, and the 
velocity distribution is given by the law of the wake 
instead of the cosine law velocity distribution. Keeping 


these differences in mind, the mean-squared error 
formula for the momentum thickness becomes 


ip 


A0„ 


= E | A0 1 
= 6 2 


(A-8) 




Substituting the law-of-the-wake velocity distribution 
(eq. (16)) into equation (A-8) and integrating gives the 
mean-square error in momentum thickness: 
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The corresponding mean-square error in the viscous 
forebody drag coefficient estimate is 




r 2 H i 2 
X\H + (2 + //)P|J 


7 2 
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(A- 10) 
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Estimating the Forebodv Drag Coefficient 
Errors 


For each configuration, the forebody pressure drag 
coefficient errors are approximated using the pressure 
coefficient standard deviations from the ten individual 
trials. For each pressure port location, 0^ , the mean and 
variance in pressure coefficients are computed as 
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(A- 13) 


S c (9) 


(A- II) 


I [C (0)-M e l 


2 I = I 


^ trials ~ ^ 


(A- 12) 


Based on equation (30), which is a pressure integral 
weighted by the cosine of the local incidence angle, 
mean-square error in forebody drag coefficient is 
computed as the sum-square of individual forebody 
pressure coefficient errors, weighted by the cosine of the 
local incidence angle. 


= X [°G COS0 1 


Estimating the Base Drag Coefficient Errors 

The base drag coefficient errors are computed in a 
similar manner as the forebody drag coefficient errors. 
However, instead of weighting the C p standard 
deviations using cos[0], the weighting procedure 
follows the scheme used to establish the base area 
curve fits (that is, the C p variances in the two ports 
along the sides of the model are weighted one-half as 
much as the C p variances in the three base-area ports): 


2 2 3 


27 

American Institute of Aeronautics and Astronautics 


References 

] Wong, Thomas J., Charles A. Hermach, John O. 
Reller, Jr., and Bruce E. Tinling, “Preliminary Studies of 
Manned Satellites — Wingless Configurations: Lifting 
Body,” NACA Conference on High-Speed 
Aerodynamics: A Compilation of the Papers Presented , 
NASA TM-X-67369, 1958, pp. 35-44. 

2 Hoerner, Sighard F., Fluid- Dynamic Drag: Practical 
Information on Aerodynamic Drag and Hydrodynamic 
Resistance , Self-published work, Library of Congress 
Card Number 64-19666, Washington, D.C., 1965. 

3 Tanner, M., “Theories for Base Pressure in 
Incompressible Steady Base Flow,” Progress in 
Aerospace Sciences , vol. 34, 1998, pp. 423-480. 

4 Rathakrishnan, E., “Effect of Splitter Plate on Bluff 
Body Drag,” AlAA Journal , vol. 37, no. 9, Sept. 1999, 
pp. 1125-1126. 

5 Saltzman, Edwin J., K. Charles Wang, and Kenneth 
W. Iliff, “Flight-Determined Subsonic Lift and Drag 
Characteristics of Seven Lifting-Body and Wing-Body 
Reentry Vehicle Configurations With Truncated Bases,” 
AIAA-99-0383, Jan. 1999. 

^Whitmore, Stephen A. and Timothy R. Moes, “A 
Base Drag Reduction Experiment on the X-33 Linear 
Aerospike SR-71 Experiment (LASRE) Flight 
Program,” AIAA-99-0277, Jan. 1999. 

7 Mills, Anthony F., Heat and Mass Transfer , Irwin 
Publishing Co., Burr Ridge, IL, 1995. 

^Adaptive Research, CFD 2000: Computational 
Fluid Dynamics System Version 2.2 User s Manual , 
Pacific-Sierra Research Corporation, 1995. 


9 Rade, Lennart and Bertil Westergren, Beta 
Mathematics Handbook: Concepts, Theorems, Methods, 
Algorithms, Formulas, Graphs, Tables , 2nd ed., CRC 
Press, Boca Raton, FL, 1990. 

10 White, Frank M. Viscous Fluid Flow , 2nd ed., 
McGraw-Hill, New York, 1991. 

n Coles, Donald, “The Law of the Wake in the 
Turbulent Boundary Layer,” Journal of Fluid 
Mechanics, vol. 1, 1955, pp. 191-226. 

12 Whitmore, Stephen A., Marco Hurtado, Jose 
Rivera, and Jonathan W. Naughton, “A Real-Time 
Method for Estimating Viscous Forebody Drag 
Coefficients,” AIAA-2000-0781, Jan. 2000 (also 
available as NASA TM-2000-2090 1 5). 

l3 Clauser, Francis H., “Turbulent Boundary Layers in 
Adverse Pressure Gradients,” Journal of Aeronautical 
Sciences , vol. 2 1 , no. 2, Jan. 1954, pp. 91-108. 

14 Krishnan, V., P. R. Viswanath, and S. Rudrakumar, 
“Effects of Riblets on Axisymmetric Base Pressure,” 
Journal of Spacecraft and Rockets , vol. 34, no. 2, 
March-April 1997, pp. 256-258. 

15 Walsh, Michael J., “Riblets as a Viscous Drag 
Reduction Technique,” AIAA Journal , vol. 21, no. 4, 
Apr. 1983, pp. 485-486. 

^Freiberger, W. F., ed.. The International Dictionary 
of Applied Mathematics , Van Nostrand Company, Inc., 
Princeton, NJ, 1960. 


28 

American Institute of Aeronautics and Astronautics 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMBNo. 0704-0188 


Puttie reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and 
maintaining the data needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this collection ol information, 
including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports. 1215 Jefferson Davis Highway, Suite 1204, Arlington, 
VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project {0704-01 B8), Washington, DC 20503 

1. AGENCY USE ONLY (Leave blank) I 2. REPORT DATE I 3. REPORT TYPE AND DATES COVERED ~ ' 


January 2001 | Technical Memorandum 


4. TITLE AND SUBTITLE 

Wind-Tunnel Investigations of Blunt-Body Drag Reduction Using 
Forebody Surface Roughness 

5. FUNDING NUMBERS 

6. AUTHOR (S) 

Stephen A. Whitmore, Stephanie Sprague, and Jonathan W. Naughton 

7 1 8-20-00-E8-53-00-b52 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSES) 

NASA Dryden Flight Research Center 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

PO. Box 273 

Edwards, California 93523-0273 

H-2439 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

Washington, DC 20546-0001 

NAS A/TM-200 1-210390 


1 1 . SUPPLEMENTARY NOTES 


Presented at 39th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, January 8-11, 2001, 
AIAA-2001-0252. Stephanie Sprague, University of Kansas. Jonathan W. Naughton, University of Wyoming. 

12a. DISTRIBUTION/ AVAILABILITY STATEMENT I 12b. DISTRIBUTION CODE 

Unclassified — Unlimited 
Subject Category 05 

This report is available at http://www.dfrc.nasa.gov/DTRS/ 

13. ABSTRACT (Maximum 200 words) 

This paper presents results of wind-tunnel tests that demonstrate a novel drag reduction technique for blunt-based 
vehicles. For these tests, the forebody roughness of a blunt-based model was modified using micomachined surface 
overlays. As forebody roughness increases, boundary layer at the model aft thickens and reduces the shearing effect of 
external flow on the separated flow behind the base region, resulting in reduced base drag. For vehicle configurations with 
large base drag, existing data predict that a small increment in forebody friction drag will result in a relatively large decrease 
in base drag. If the added increment in forebody skin drag is optimized with respect to base drag, reducing the total drag 
of the configuration is possible. The wind-tunnel tests results conclusively demonstrate the existence of a forebody drag- 
base drag optimal point. The data demonstrate that the base drag coefficient corresponding to the drag minimum lies 
between 0.225 and 0.275, referenced to the base area. Most importantly, the data show a drag reduction of approximately 
15 percent when the drag optimum is reached. When this drag reduction is scaled to the X-33 base area, drag savings 
approaching 45,000 N (10,000 lbf) can be realized. 


14. SUBJECT TERMS 15. NUMBER OF PAGES 

Base drag. Drag reduction. Reusable launch vehicle, Skin friction. Wind tunnel ^4 

16. PRICE CODE 

A03 

17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION OF ABSTRACT 

OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified Unlimited 

Standard Form 298 (Rev, 2-89) 

Prescribed by ANSI Sid Z39-1B 
298-102 


NSN 7540-01-280-5500 























National Aeronautics and 
Space Administration 
Code JTT 

Washington, D.C. 20546-0001 
USA 

Official Business 
Penalty for Private Use, $300 


SPECIAL FOURTH-CLASS RATE 
POSTAGE AND FEES PAID 

MAQA 

PERMIT No G27 




'k 


POSTMASTER: If Undeliverable (Section 158 
Postal manual) Do Not Return 



